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ABSTRACT 

We report on a deep search for radio pulsations toward five unidentified ASCA X-ray sources coincident 
with EGRET 7-ray sources. This search has led to the discovery of a young and energetic pulsar using 
data obtained with the new Wideband Arecibo Pulsar Processor. PSR J2021+365I is likely associated 
with the X-ray source AX J2021. 1+3651, which in turn is likely associated with the COS B high energy 
7-ray source 2CG 075+00, also known as GeV J2020+3658 or 3EG J2021+37I6. PSR J2021+3651 has 
a rotation period P = 104 ms and P = 9.6 x 10 -14 , implying a characteristic age r c ~ 17 kyr and a 
spin-down luminosity E ~ 3.4 x 10 36 ergs s _1 . The dispersion measure DM~ 371 pc cm' 3 is by far 
the highest of any observed pulsar in the Galactic longitude range 55° < I < 80° . This DM suggests a 
distance d > 10 kpc, and a high 7-ray efficiency of ~ 15%, but the true distance may be closer if there is 
a significant contribution to the DM from excess gas in the Cygnus region. The implied luminosity of the 
associated X-ray source suggests the X-ray emission is dominated by a pulsar wind nebula unresolved 
by ASCA. 

Subject headings: pulsars: general — pulsars: individual (PSR J2021+3651) — stars: neutron — 

X-rays: individual (AX J2021. 1+3651) — gamma-rays: individual (GeV 2020+3658) 



1. INTRODUCTION 

The majority of high energy 7-ray sources observed by 
EGRET and other telescopes have long escaped identifica- 
tion with lower energy counterparts (Hartman et al. 1999). 
Young pulsars remain the only Galactic source class (other 
than the Sun) unambiguously shown to emit radiation in 
the 100 MeV - 10 GeV range (Thompson 2001). It is 
likely that many of the unidentified 7-ray sources at low 
Galactic latitudes are young pulsars as well. Many of these 
sources have characteristics similar to those of the known 
7-ray pulsars, but have no known pulsars within their er- 
ror boxes. This fact, along with modelling of the multi- 
wavelength pulse profiles and the still singular example of 
Geminga (Halpern & Holt 1992; Helfand 1994), has led 
to the suggestion that a large fraction of the radio beams 
from 7-ray sources will miss the Earth and appear radio 
quiet (Romani 1996). 

Recently, a number of young pulsars coincident with 
known 7-ray sources have been discovered (D Amico et al. 
2001; Camilo et al. 2001). These new discoveries are 
largely a result of greater sensitivity to pulsars with high 
dispersion measures (DM) obtainable with newer pulsar 
backends such as the Parkes multibeam system (Manch- 
ester et al. 2001). The recent detection of a young radio 



pulsar in the supernova remnant 3C58 with a 1400 MHz 
flux density of only ~ 50 /iJy (Camilo et al. 2002) suggests 
many more faint radio pulsars await discovery in deep, tar- 
geted searches. 

A major stumbling block in the identification of the 
EGRET sources is their large positional uncertainty, which 
can be greater than 1° across. We approach this problem 
by targeting potential hard X-ray counterparts, whose size 
and positional uncertainty are much smaller than the typ- 
ical single dish radio beam. Using as our guide the ASCA 
catalog of potential X-ray counterparts of GeV sources 
(based on the Lamb & Macomb (1997) catalog of sources 
with significant flux above 1 GeV) by Roberts, Romani, 
& Kawai (2001) (hereafter, RRK), we have searched five 
X-ray sources for radio pulsations using the 305-m Arecibo 
telescope and the 64-m Parkes telescope (see Table 1). 
Previous searches of these targets were limited. In par- 
ticular, two of the three sources observed at Parkes (AX 
J1418. 7-6058 and AX J1809.8-2332) were not previously 
the subject of any directed search and were only observed 
as a matter of course during the Parkes Multibeam Galac- 
tic Plane Survey (Manchester et al. 2001). A survey of 
EGRET sources by Nice & Sayer (1997) looked at two 
of the sources searched here (AX J1826. 1 — 1300 and AX 
J2021. 1+3651) with a limiting flux density for slow pulsars 
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of 0.5 - 1.0 mJy at frequencies of 370 and 1390 MHz, but 
found no new pulsars. Our search has led to the discovery 
of one young and energetic pulsar, PSR J2021+3658. We 
argue that it is a likely counterpart to AX J2021. 1+3651 
and GeV J2020+3658 / 2CG 075+00. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Observations 

On 2002 January 30 and 31, we observed the only 
two unidentified sources in the RRK catalog visible from 
the Arccibo radio telescope, AX J1907. 4+0549 and AX 
J2021. 1+3651, using the Wideband Arecibo Pulsar Pro- 
cessor (WAPP). The WAPP is a fast-dump digital cor- 
relator with adjustable bandwidth (50 or 100 MHz) and 
variable numbers of lags and sample times (for details see 
Dowd, Sisk, & Hagen 2000). Our observations were made 
at 1.4 GHz with 100 MHz of bandwidth and summed po- 
larizations. The observational parameters are summarized 
in Table 1. The 16-bit samples were written to a disk array 
and then tranfered to magnetic tape for later analysis. 

On 2001 February 11 — 15, the three extended hard X- 
ray sources listed by RRK as potential pulsar wind neb- 
ulae, AX J1418. 7-6058 (the Rabbit), AX J1809.8-2333, 
and AX J1826. 1 — 1300, were searched for radio pulsations 
with the Multibeam receiver on the Parkes radio telescope. 
Each source was observed once at a central observing fre- 
quency of 1390 MHz with a 512 channel filterbank spec- 
trometer covering 256 MHz of bandwith, and once at a 
central observing frequency of 1373 MHz with 96 chan- 
nels and 288 MHz of bandwith (see Table 1). During each 
observation, signals from each channel were square-law de- 
tected and added in polarization pairs before undergoing 
high-pass filtering. The signals were one-bit digitized ev- 
ery 0.25 ms and recorded onto magnetic tape for later 
analysis. 

2.2. Analysis 

Analysis of Arecibo observations was done using the 
PRESTO software suite (Ransom 2001) by first removing 
obvious narrow band and/or short duration interference 
in both the time and frequency domains. We then dedis- 
persed the data at 500 trial DMs between 10 and 510 pc 
cm- 3 for AX J2021. 1+3651 and 540 trial DMs between 
and 2695 pc cm" 3 for AX J1907.4+0549. Employing 
harmonic summing, the FFTs of each time series were 
searched, and interesting candidates were folded over a 
fine grid in DM, period, and period derivative space to 
optimize the signal-to-noise. 

The Parkes observations were analyzed using standard 
pulse search software and a similar procedure by searching 
the 96 channel data at 279 trial DMs ranging from to 
1477 pc cm' 3 and the 512 channel data at 501 trial DMs 
ranging from to 670 pc cm~ 3 . We tested the system 
by observing a known bright pulsar (PSR B1124— 60) for 
300 s, which was clearly detected in the processing. Re- 
analysis of the data using PRESTO has not revealed any 
new candidates. 

3. RESULTS 

3.1. PSR 32021+3651 

2 See http://pulsar.princeton.edu/tempo 



A new highly dispersed 104-ms pulsar was detected in 
the Arecibo observations made of AX J2021. 1+3651; it is 
clearly visible in both of the original search observations, 
and represents a ~ 40 a detection in the longest observa- 
tion. The pulse profile is shown in Figure 1. 

A subsequent series of 7 observations performed between 
MJD 52405—52416 allowed us to determine a phase con- 
nected solution for some of the pulsar parameters. In- 
tegrated pulse profiles from these observations were con- 
volved with a template profile to extract 12 topocen- 
tric times of arrival (TOA). Using TEMPO 2 and adopting 
the ROSAT position for the pulsar (§4.1), the topocen- 
tric TOAs were converted to TOAs at the solar system 
barycenter at infinite frequency and fit simultaneously for 
pulsar period, period derivative, and DM, with a residual 
rms of 91 fis. The measured and derived parameters for 
this pulsar are listed in Table 2. 

3.2. Non-detections 

No convincing pulsar candidates were detected in any 
of the search observations conducted at Parkes. We esti- 
mate upper limits of S = 0.08 mJy at 1.4 GHz for pulse 
periods P > 10 ms for most of the 512 channel observa- 
tions. A comparable sensitivity was obtained for long pe- 
riods in the 96 channel observations. For DMs larger than 
about 100 pc cm~ 3 , the sensitivity to fast pulsars (P < 
50 ms) is significantly degraded in the 96 channel system. 
These sensitivity limits were estimated using a sensitiv- 
ity modeling technique described in detail elsewhere (e.g., 
Manchester et al. 2001). 

Likewise, extensive searching of AX J1907. 4+0549 
yielded no convincing pulsar candidates. We estimate an 
upper limit of S ~ 0.02 mJy at 1.4 GHz for a long period 
pulsar assuming a 10% duty cycle. 

4. DISCUSSION 

4.1. PSR 32021+3651 

Our search targeted AX J2021. 1+3651, which was iden- 
tified as a potential high-energy counterpart to GeV 
J2020+3658 by RRK. The X-ray source is near the ASCA 
field edge and so the positional uncertainty from ASCA is 
> 1' (Gotthelf et al. 2000). A subsequent search of the 
R OS A T All-Sky Survey Faint Source catalog (Voges et al. 
2000) revealed the source 1RXS J202104.5+365127 with a 
smaller positional error of 24". Given the rarity of such 
young, energetic pulsars and the small size of the Arecibo 
beam (3' at FWHM), an association with the X-ray source 
is highly probable. 

The DM of PSR J2021+3651 is by far the highest known 
in the Galactic longitude range 55° < I < 80° which is 
mainly an inter-spiral arm direction. The Taylor & Cordes 
(1993) model gives a distance of ~ 19 kpc, well beyond the 
last spiral arm used in the model. A revised model is cur- 
rently in preparation by Cordes and Lazio which includes 
an outer spiral arm at d ~ 10 kpc. Placing the pulsar at 
the far edge of this outer arm still does not account for all 
the observed dispersion; and it is possible that there are 
further contributions from clouds in the Cygnus region, 
where there is known to be excess gas at d ~ 1.5 kpc (J. 
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Cordes, private communication). However, there are no 
obvious HII regions within the Arccibo beam seen in ei- 
ther Very Large Array (VLA) 20-cm radio or Midcourse 
Space Experiment (MSX) 8.3 /im images (available from 
the NASA/IPAC Infrared Science Archive). 

The high DM is somewhat surprising given the X-ray 
absorption quoted by RRK, n H =(5.0±0.25) x 10 21 cm~ 2 , 
where the errors represent the 90% confidence region. The 
total Galactic HI column density in this direction as es- 
timated from the FTOOL nh, which uses the HI map of 
Dickey & Lockman (1990), is 1.2 x 10 22 cm" 2 . This should 
be a good approximation if the source is truly at the far 
edge of the outer spiral arm. Noting that the ASCA im- 
age shows faint, softer emission in the region (Figure 2), 
and given the likely possibility of either associated thermal 
X-ray flux from a supernova remnant or a nearby massive 
star, we fit the ASCA spectrum of RRK, adding a thermal 
component to the absorbed power-law model. Account- 
ing for ~ 4% of the photon flux with a MEKAL thermal 
plasma model of temperature kT <~ 0.1 keV in XSPEC 
(Arnaud 1996) statistically improves the fit (F-test chance 
probability of 2.5%). The best-fit absorption for this three 
component model is njj=7.6xl0 21 cm~ 2 with a 90% confi- 
dence region of (4.1 - 12.3) x 10 21 cm~ 2 , consistent with the 
total Galactic column density. The best-fit photon index 
is T = 1.86, still consistent with the 1.47 - 2.01 range in 
RRK derived from the simple absorbed power-law model. 
Hence the X-ray absorption does not force us to adopt a 
smaller distance than is suggested by the DM. 

For a distance d w = rf/10 kpc, the inferred isotropic X- 
ray luminosity L x = 4.8 x 10 34 (i 2 (2 - 10 keV). The X-ray 
efficiency rjx = Lx/E is 0.01g? 2 . Compared to the total 
pulsar plus nebula X-ray luminosity of other spin-powered 
pulsars this is somewhat high, but within the observed 
scatter (Possenti et al. 2002; Chevalier 2000). 

The pulsar's positional coincidence with the error box 
of the hard spectrum, low variability EGRET 7-ray source 
GeV J2020+3658 coupled with the high inferred spin-down 
luminosity strongly suggests this pulsar emits pulsed 7- 
rays. Unfortunately, confirming this by folding archival 
EGRET data is problematic due to the likelihood of sig- 
nificant past timing noise and glitches, which make the 
back-extrapolation of the rotational ephemeris uncertain. 
RRK noted that the chance probability of an X-ray source 
as bright as AX J2021. 1+3651 in the EGRET error box 
was ~ 10%, but the nearby Wolf-Rayet star WR141 was 
equally bright in X-rays and also a potential 7-ray emit- 
ter. However, young pulsars remain the only firmly estab- 
lished class of Galactic EGRET sources. The known 7-ray 
pulsars cluster at the top of pulsar lists rank-ordered by 
spin-down flux E/d 2 , with 7-ray efficiencies ry 7 = L 1 /E 
mostly between 0.001 and 0.03 (assuming 1 sr beaming) 
with a tendency to increase with pulsar age (Thompson 
et al. 1999). The exception is PSR B1055-52, with an 
apparent 7-ray efficiency 7y 7 ~ 0.2 given its nominal DM 
distance of 1.5 kpc. The inferred 7-ray efficiency for PSR 
J2021+3651 is r/ 7 = 0.15rf 2 in the 100 MeV to 10 GeV 



range. If the pulsar is located within the Perseus arm 
at a distance of 5 kpc, then the inferred X-ray and 7- 
ray luminosities would be fairly typical of the other pul- 
sars with Vela-like spin-down luminosities. While there 
is currently no observational evidence for a distance this 
close, increased DM from an intervening source in this rel- 
atively crowded direction would not be surprising. We 
note that the DM derived distance for another young pul- 
sar recently discovered within an EGRET error box, PSR 
J2229+6114, also leads to an anomalously high inferred 
7-ray efficiency (Halpern et al. 2001). 

4.2. Upper Limits Toward the Other Sources 

Determining the fraction of radio-quiet versus radio- 
loud pulsars is important for our understanding of 7-ray 
pulsar emission mechanisms. The two leading classes of 
emission models, the outer-gap (Romani 1996) and polar- 
cap (Daughcrty & Harding 1996) models, make very dif- 
ferent estimates of the fraction of 7-ray pulsars that should 
be seen at radio energies. Out of the 25 brightest sources 
above 1 GeV not associated with blazars, <~ 10 are now 
known to either be energetic radio pulsars or contain such 
pulsars within their error boxes. Searching the brightest 
unidentified X-ray sources in five GeV error boxes, we de- 
tected radio pulsations at the ~ 0.1 mJy level (similar to 
the limiting sensitivity of the Parkes observations) from 
one of these with Arecibo. This is well below the average 
flux level expected for typical radio luminosities of young 
pulsars (Brazier & Johnston 1999) and distances to star 
forming regions statistically associated with 7-ray sources 
(Yadigaroglu & Romani 1997). Two of the sources ob- 
served with Parkes, AX J1418.7-6058 (the Rabbit) and 
AX J1809.8— 2333, have radio and X-ray properties that 
clearly identify them as pulsar wind nebulae (Roberts et al. 
1999; Braje et al. 2002), and the third, AX J1826. 1-1300, 
is an extended hard X-ray source that has few other source 
class options. Therefore, all three remain viable candi- 
dates for 7-ray loud, radio-quiet pulsars. Out of this same 
sample of 25 bright GeV sources, the total number of 
reasonable candidate neutron stars within the 7-ray er- 
ror boxes which have now been searched deeply for radio- 
pulsations without success is ~ 7. A current "best guess" 
fraction of radio-loud 7-ray pulsars of <~ 1/2 falls in be- 
tween the predictions of the two main competing models. 
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Table 1 



Observational Parameters 



Source 


RA (J2000) 


DEC (.72000) 


Epoch 


Telescope 


fc 8 

(MHz) 




N ch c 


t d 


Tint 6 


S ml „ f 


Name 


hh:mm:ss.s 


dd:mm:ss 


(MJD) 




(MHz) 




(pa) 


(s) 


(mJy) 


AX J1418. 7-6058 / GeV J1417-6100 


14:18:41.5 


-60:58:11 


51951.63 


Parkcs 


1390 


256 


512 


250 


16900 


0.08 






51953.67 


Parkcs 


1373 


288 


96 


250 


16900 


0.08 


AX J1809.8-2332 / GeV J1809-2327 


18:09:50.2 


-23:32:23 


51951.83 


Parkcs 


1390 


256 


512 


250 


16900 


0.08 






51954.77 


Parkcs 


1373 


288 


96 


250 


16900 


0.08 


AX J1826.1-1300 / GeV J1825-1310 


18:26:04.9 


-12:59:48 


51952.83 


Parkcs 


1390 


256 


512 


250 


13234 


0.09 






51955.83 


Parkcs 


1373 


288 


96 


250 


14832 


0.09 


AX .71907. 4+0549 / GeV .71907+0557 


19:07:21.3 


+05:49:14 


52305.58 


Arccibo 


1425 


100 


512 


200 


6480 


0.02 


AX J2021. 1+3651 / GeV J2020+3658 


20:21:07.8 


+36:51:19 


52304.67 


Arccibo 


1425 


100 


512 


200 


1627 


0.04 






52305.66 


Arccibo 


1425 


100 


256 


200 


3000 


0.03 



a Ccntral observing frequency. 

b Total bandwidth of the observation. 

c Numbcr of frequency channels. 

d Sampling time. 

°Lcngth of the observation. 

f Flux density sensitivity limit. 
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Table 2 

Measured and Derived Parameters for PSR J2021+3651 



Parameter 


Value 


Right ascension a (J2000) 


20 h 21 m 04 s . 5 a 


Declination S (J2000) 


+36° 51' 27". a 


Galactic latitude I 


75.° 23 


Galactic longitude b 


+0.°11 


Pulse period P (s) 


0.10372222480(17) 


Period derivative P 


9.563(48) x 10" 14 


Pulse frequency v (s^ 1 ) 


9.641135271(16) 


Frequency derivative v (s~ 2 ) 


-8.889(45) x 10~ 12 


Epoch (MJD) 


52407.389 


Dispersion Measure DM (pc cm~ 3 ) 


371(3) 


Pulse width at 50% of peak W50 (ms) 


9.9 


Pulse width at 10% of peak wio (ms) 


18 


Flux density at 1425 MHz (mjy) 


~ 0.1 


Spin-down luminosity P b (erg s _1 ) 


3.4 x 10 36 


Surface dipole magnetic field B C (G) 


3.2 x 10 12 


Characteristic Age r c = ^P/P (kyr) 


17 



Note. — Figures in parentheses represent uncertainty in the 
least-significant digits quoted, equal to 3 times errors given by 
TEMPO. 

a Coordinates of the ROSAT X-ray source 1RXS 
J202104. 5+365127 with estimated positional error of 24". 

b E = ±tt 2 IP/P 3 with / = 10 45 g cm 2 . 

c Assuming standard magnetic dipole spindown: B — 3.2 x 
10 19 (PP) 1/2 Gauss (Manchester fc Taylor 1977). 
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Fig. 1 — 



Pulse Phase 

1.4 GHz pulse profile for PSR J2021+3651 from the MJD 52305 observation. The error bar represents the 1 <r uncertainty. 
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Fig. 2. — ASCA GIS 2-10 keV image of the 7-ray source region. The contours are the 68%, 95%, and 99% confidence regions of the 7-ray 
source position, derived from the > 1 GeV photons (Roberts, Romani, & Kawai 2001). WR 141 is a Wolf-Rayet star also in the field of view. 
The circle centered on AX J2021. 1+3651 indicates the size of the 3' (FWHM) Arecibo beam. 



